DC four point resistance of a double barrier quantum pump 
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We investigate the behavior of the dc voltage drop in a periodically driven double barrier structure 
(DBS) sensed by voltages probes that are weakly coupled to the system. We find that the four 
terminal resistance Ru measured with the probes located outside the DBS results identical to the 
resistance measured in the same structure under a stationary bias voltage difference between left 
and right reservoirs. This result, valid beyond the adiabatic pumping regime, can be taken as an 
indication of the universal character of Ru as a measure of the resistive properties of a sample, 
irrespectively of the mechanism used to induce the transport. 
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Quantum transport induced by time dependent fields 
attracts presently an impressive amount of research. A 
phase coherent conductor subjected to two periodically 
varying voltages becomes a paradigmatic example of a 
quantum pump, in which a dc current can be generated 
in the absence of a net external hias^^^^^. 

After the works of Landauer and Buttiker,.^'- the four 
point resistance Ra is considered as the proper measure 
of the genuine resistive behavior of a mesoscopic sample, 
free from the effects of the contact resistance. Several 
theoretical works have been devoted to study the de- 
tails of the voltage drop between the contacts in systems 
where the transport is induced by means of a stationary 
dc voltage bias^. Furthermore, the striking feature that 
i?4t can be negative in a coherent conductor has been ex- 
perimentally observed in semiconductors^ and in carbon 
nanotubes^. However, the behaviour of R^t in the case 
of a quantum pump has not been so far analyzed. 

The aim of the present work is to investigate to what 
an extent the concept of Rn could depend on the un- 
derlying driving mechanism. To this end, we consider 
as a model of the quantum pump, a quantum wire cou- 
pled to left and right reservoirs at a fixed chemical po- 
tential and with two narrow gates to which oscillating 
voltages are applied with a phase-lag. This set up (see 
FiglT] upper plot) mimics the actual double barrier struc- 
ture (DBS) used in RefQ where two of such ac potentials 
were applied at the walls confining a quantum dot. Ex- 
perimentally, the dc response is actually inferred from 
the measurement of the voltage drop between two ex- 
tra probes, one located at the left and the other at the 
right of the DBS. Accordingly, we consider non- invasive 
voltage probes weakly coupled to the wire. The chemi- 
cal potential fii {i — P, P ) of each probe is adjusted to 
maintain zero net current through the respective contact. 
Our goal is to investigate the dc four terminal resistance 
for the quantum pump defined as the dc voltage drop 
Afippi sensed by the two probes P, P' connected at two 
arbitrary points along the sample, divided by the dc com- 
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FIG. 1: (Color on-line) Scheme of the set up. The central 
device is a wire with two barriers of height Eb connected to 
L and R reservoirs. Two voltage probes P and P sense the 
voltage drop. Upper plot: Pumping set up in which a dc 
current J'*'^ is induced by two local ac voltages. The L and R 
reservoirs are at the same chemical potential fj,. Lower plot: 
Stationary set up, in which the current J" is induced by a dc 
voltage difference 14 between L and R reservoirs. See text for 
more details. 



ponent of the current J'^'^ flowing through the device : 



We will also compare this quantity with the four terminal 
resistance obtained when the current through the DBS 
is induced by a slight stationary voltage difference Vs 
between the left (L) and right (R) reservoirs, as indicated 
in the lower panel of Fig[TJ 

For sake of clarity we start considering just one voltage 
probe, which is modeled as a third reservoir coupled to 
the central system at the position P (the extension to 
more probes is trivial). The corresponding Hamiltonian 
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for the full system reads: 

H ^ Hieads+Hp + Hc{t)-WL{a{ci+H.c.) 

—WRia^RCN + H.c.) - wp{a)pCp + H.c.) (2) 

with Hc{t) denoting the Hamiltonian for the central 
piece that we model as a Id tight-binding chain of length 
iV with two barriers located at sites A and B: 

N 

1=1 

N 

-Wh ^(c|c/+i + H.c) + V cos.{^ot)c^gCB, (3) 
1=1 

with Wh the hopping parameter and the profile ea — 
Eb — Eb, El — 0,1 — 1, . . . , N A, B, defining a double 
barrier structure. The time dependent ac potentials act 
locally at the position of the barriers and have amplitude 
V, frequency JIq oscillate with a phase difference S. 
We denote with Hi^ads the Hamiltonians of two semi- 
infinite tight-binding chains with hopping wi at the same 
chemical potential /x, which play the role of the L and 
R reservoirs. These two leads are connected to the cen- 
tral device at sites 1, N respectively. Similarly, Hp is the 
hamiltonian of the voltage probe P that we also model 
as a particle reservoir with a chemical potential /ip that 
is fixed to satisfy the condition of net zero dc current 
through its contact to the central system^. The con- 
tacts between the central system and the L and R leads 
and the probe P are described by the last three terms of 
Eq.(I21), where the fermionic operators Oq (a = L,R,P) 
denote degrees of freedom for the L, R and P reservoirs 
respectively. 

We employ the formalism of Keldysh non equilibrium 
Green's functions, which is a convenient tool in trans- 
port theory on multiterminal structures driven by time- 
periodic fields. Following Ref 4 we use the Floquet rep- 
resentation G^'i,{t,uj) = J2^^^^e~'''^°^gi^i,{k,uj), where 
Gj^i, {t, oj) is the Fourier transform with respect to t — t' 
of the retarded Green's function. The dc component of 
the charge current flowing through the contact between 
the central system and the probe P, can be written (in 
units of e/h) as^: 

Jp= E / l^{raHTp{u; + kno) 

a=L,P,Rk=^oo-'-°° 

\gi,^i^ {k, u;)\^ iUu;) - fp{u + mo)]}, (4) 



where la are the sites of the central system at which 
the reservoirs a = L,R,P are attached, while Ta{uj) = 
\wa\'^Pa{i^) is the spectral function associated to the self- 
energies due to the coupling to these reservoirs, Pq(ijj) 
is the corresponding density of states and — 



1/0 



-|- 1) is the Fermi function of the reservoir 



that satisfies Jp'^ = in the above expression. It is, how- 
ever, instructive to analyze first the case of low driving 
frequency fio and small pumping amplitude V, which cor- 
responds to the so called adiabatic pumping regime^. As 
we already mentioned, we are considering "non-invasive 
probes" . This corresponds to probes weakly coupled to 
the central system, in such a way that they do not in- 
troduce neither inelastic nor elastic scattering processes 
for the electronic propagation between L and R reser- 
voirs. Below we derive an analytical expression, valid 
under these conditions, for the voltage profile /ip. 

For low V, the perturbative solution of the Dyson's 
equation up to the second order in this parameter con- 
tains only the following Floquet components that con- 
tribute to the dc current^: 

+e±'^G°s(c.Tr!o)G^,i'M]. (5) 

For weak coupling to the probes, Eq.Q is evaluated 
with Green's functions up to the 1st order in wp. This 
corresponds to consider the functions GW,{uj) in as 
the equilibrium retarded Green's functions of the cen- 
tral system attached only to the L and R reservoirs. 
For perfect matching to the reservoirs (wp = wp — 
wi — Wh) and for barriers with low amplitude Ep < Wh, 
these functions can be written in the following simple 
form: G?,, (c.) = 5^,^^) + ^^5 E,=A.s with 

gii,{0) = ie-'\^-^'\y{2whSm0), being lu = 2u;,, cos 6'^^. 
Using them to evaluate ([5]) , substituting the result in Q 
and considering the adiabatic (oc flo) contribution in the 
resulting Jp*^ we get two different results depending on 
the place at which the probe is connected: 

/ip = fi±noV'^ sinS[a°{kp) 

+i?B/3°(fcp)], Xp>XB,Xp<XA, 

lip = fi + rioV'^ sinSla^kp, xp) 

+EBl3\kp,xp)], XA<xp<XB, (6) 

with Xj{j = A, B, P) denoting the position of the barriers 
and probe in units of the lattice parameter of the tight- 
binding model. The upper and lower sign of the first 
identity corresponds, respectively, to the voltage probe 
located at the left {xp < xa) and right {xp > xp) side 
of the DBS, while the second identity corresponds to the 
voltage probe located between the two barriers. We have 
defined the Fermi vector (in units of the lattice parame- 
ter) as kp = 6{fi) as well as the following functions: 



a°(fcp) — 



a, which we assume to be at the temperature = 0. 

The voltage profile sensed by the probe can be exactly 
evaluated under general conditions from the solution fip 



P°{kF) 

a^{kF, Xp) 
f3\kF,xp) 



sin[2kp{xA — xp)] 

2{wh sin fcp)^ 
sin^[fcp(xA - Xp)] 



A{whSaikp)'^ ' 
= sin[fcp(2a;p — xa — XB)]a°{kp), 
^ii\n[kp{2xp~XA~XB)]l3°{kp), (7) 

where the superscripts o, («) stress that the probe senses 
points outside (inside) the region where all the scattering 
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FIG. 2: (Color on-line) Local voltage fip sensed by the volt- 
age probe P as a function of the probe position xp along a 
DBS composed by A'^ = 50 sites with two barriers of height 
Eb = 0.2 located at xa ~ 5 and a;_g = 45 as indicated by the 
vertical dashed lines. The pumping parameters are V — 0.01, 
Slo = 0.01 and S = 7r/2. Red squares correspond to Eq.®, 
the analytical solution for the adiabatic pumping regime and 
a weakly connected probe. Black circles correspond to the 
exact numerical solution obtained equating Eq. Q to zero 
with wp = 0.01. The chemical potential is = 0.2, which 
corresponds to kp — 1.47. 



processes (dynamical as well as stationary) take place. 
In the simple model we are considering, with a perfect 
matching between the central system and the reservoirs, 
this coincides with the DBS, as indicated in Fig. 1. Equa- 
tions ^ and ([7]) tell us that the local voltage sensed 
by a probe is constant outside the region defined by 
the DBS, while it presents the characteristic pattern of 
Friedel oscillations^ with a period 2kp at positions lying 
between the two oscillating barriers. Fig. [5] shows the 
benchmark of the analytical result, Eq.(IS]), against the 
exact voltage profile obtained numerically from Eq. (|10p 
in the regime of weak V, and wp, and a moderate Ep- 
A good agreement of the qualitative behavior is observed. 
In particular, the exact profile fip exhibits Friedel oscilla- 
tions with period 2kp as a function of the probe position 
Xp as predicted by Eq.([S]) and only a slight disagreement 
is found in the amplitude of the envelope function. 

To the lowest order of perturbation in the coupling 
constant wp, the effect of an additional second voltage 
probe P' is completely uncorrelated from the first one, 
since the associated interference effects involve second 
order processes in wp. At this level of approximation 
let us call upr the local voltage sensed by the additional 
probe at P' and Afippi = fip> — fip the corresponding 
voltage drop. In a set up with the probe P (P') located at 
the left (right) side of the DBS, the voltage drop between 
both probes is from Eq.(l6]): 

AVpp- = 2r2oV^^ sin(5[a°(fcp) + EBl3°{kF)]. (8) 

Another possible measurement corresponds to locate the 
voltage probes P and P' inside the DBS. In this case, the 
voltage drop between the two probes explicitly depends 



on the probe positions xp and xpi as follows: 



a''{kF)+EB[3°{kF) 



AVpp- =2r!oF2sin5 . , 

[ s\n{kp{XA -xb)) 

cos{A;p[(a::p — xa) + [xp' ~ xb)]} sin[A;p(a;p' — xp)] 



(9) 



where, as before, we have employed the superscripts o, [i) 
to distinguish configurations with the probes outside (in- 
side) the DBS. 

Under the conditions assumed in the derivation of 
Eq.([6]), i.e. low y, rio, wp and Ep, the dc current flowing 
through the DBS reads: 



jdc ^ 2Tlr%noV^ sin (5 



a''{kF) + EBP°ikF) 



{2wh sin fcp)2 



(10) 



with = r° (/i), {a = L,R). We can now compute the 
dc four terminal resistance Rn in an adiabatic weakly 
driven pumping process. For a set up in which the probes 
are located outside the DBS {xp < xa,xb < xpi) it 
reads: 



■ ^^pp' 

Jdc 



{2wii sin k} 



(11) 



while it is: 



A>pp/ 

Jdc 

sin[fcp(a;p/ - xp)] 
SU\[kF{XA - xb)\ 
X cos[fcp(xp — XA + xp' — Xb)] , (12) 

for a set up in which the probes are located inside the 
DBS {xa < Xp < xp' < Xp)- 

We now turn to compare the value of i?4t obtained for 
the quantum pump with the resistance of the same DBS 
when the transport is induced through a stationary bias 
voltage Vs established by a difference in the electrochem- 
ical potentials of L and R reservoirs fj,L = fi + eVs/2 
and ftp = /i — eVs/2, as is depicted in the lower panel of 
FiglH Under the conditions of non-invasive probes and 
for linear response in Vs, it is possible to implement the 
same kind of perturbative procedure as before, but for 
the stationary Green's functions. Assuming again perfect 
matching between the DBS and the L and R reservoirs 
and Eb < Wh, we derive the following simple expression 
for the current flowing through the device, 



r 



{2wh sin fcp)^ 



(13) 



which is the stationary counterpart of Ea. (jlO[) . ^From 
the above expression we compute Vg/J'^ and arrive im- 
mediately to the important identity: 



^4* = J7 



(14) 



which tells that the dc four point resistance measured in 
the quantum pump when the probes are connected out- 
side the DBS i?4j (Eq. pTjl ) exactly coincides with the 
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total resistance of the structure measured under station- 
ary bias, provided that the driving condition corresponds 
to linear response in the stationary setup and the adia- 
batic regime in the pumping setup. At this point, it is 
important to recall that in the outside configuration the 
two probes enclose the whole region where all the scat- 
tering processes and, therefore, the full voltage drop Vg 
applied in the stationary setup takes place. 

On the other hand, the counterpart of Eq. ^2]) for the 
stationary configuration reads: 

=^4* ' ^^^^ 

which means that inside the DBS the dc resistance for the 
pumping setup differs from that under stationary driving 
just in a geometrical factor. 

Besides the relevance of the analytical results it is valu- 
able to analyze the response of the system beyond the 
adiabatic pumping condition. For that purpose we have 
performed extensive numerical calculations of the dc cur- 
rent flowing through the system J'^'^ and the potential 
drop sensed by the voltage probes along the DBS, as a 
function of the pumping frequency Qq. In the left panel 
of Fig.© we plot the chemical potential /ip as a func- 
tion of r^o for different locations of the probe xp. In 
addition, in the right panel we have plotted the voltage 
drop A/i = fipi — fip between two probes located outside 
the DBS (circles) and inside the DBS (squares). 

When the frequency ^Iq is close to the energy differ- 
ence between two neighbouring levels of the DBS, the 
latter become mixed by the pumping potential which 
causes an inversion in the sign of the dc current. A 
rough estimate for the frequency at which such a res- 
onant condition is achieved in the example of Fig. 3 
casts r2o ~ 0.22. In good agreement, we find an inversion 
in the sign of A°fippi for probes connected outside the 
DBS at flo ^ 0.25 (plots in circles of Fig. 3). Moreover, 
the voltage drop between two points outside the DBS in 
this case results identical to the product of the dc cur- 
rent J'^'^ times the value of the four point resistance i?!^ 
obtained in the adiabatic pumping regime in Eg. pip . In 
other words, our results indicate that even for pumping 



frequencies flo beyond the adiabatic regime, it is still pos- 
sible to unambiguously define R^^ as the value obtained 
under the adiabatic approximation fEg. ljlip ). On the 
other hand, the voltage drop measured inside the DBS 
coincides with the product of R\f times J'^'^ only within 
the adiabatic pumping regime, i.e when A''fippi depends 
linearly on Qq- 

In conclusion we have shown that the four terminal 
resistance i?4j measured in a pumping set up for probes 
located outside the DBS coincides with the total resis- 
tance of the structure measured under stationary bias. 
This result can be taken as an indication of the universal 
character of Rl^ as a concept to characterize the resis- 
tive properties of a system. Our calculation could be 
extended to the case where the probes themselves are 
subjected to ac voltages. 
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FIG. 3: (Color on-line) Left panel: Exact local voltage /ip 
sensed at sp = 2, 14 (empty symbols) and Xp/ = 26, 38 (solid 
symbols) , as a function of the pumping frequency Qo obtained 
numerically equating Eq. Q to zero. Circles (squares) corre- 
spond to points outside (inside) the DBS. Parameters are: 
fj, = 0.4, N = 40, Eb = 0.2, XA = 10 and xb = 30. 
Right panel: The corresponding potential drops A°''fippi. 
The products R'^^J'^'^ (solid and dashed lines, respectively) 
between the exact current J'*'^ and the resistances evaluated 
from Eqs. (|lip and (|12|l . respectively. The vertical dotted 
line indicates the frequency of the pumping at which current 
inversion occurs. 
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